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Abstract 

It is shown that solutions to the 2nd order BFKL eigenvalue equation exist for arbitrary 
large real values of the complex angular momentum j. This corresponds to a cut in the 
complex j plane along the whole real axis, and it makes the use of the complex angular 
momentum variable for the calculation of the high-energy behavior somewhat questionable. 
The eigenfunctions contain non-perturbative pieces which behave as exp(—l/as6) and have 
no counterpart in the leading-log BFKL equation. The high-energy behavior of the 2nd 
order BFKL Green function as found by other authors, is reproduced by excluding these 
non-perturbative pieces of the eigenfunctions. 
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1 Introduction 


Recently corrections of the 2nd order in the conpling constant to the hard pomeron 
(BFKL) equation have been calculated |I|. For the intercept they are accompanied by 
a large negative coefficient Bl], which makes its prediction reliable only for extremely 
low coupling constants and translates into extremely large relevant mass scales. Al¬ 
though this conclusion has been questioned in |^, different calculations iii show 
features of this approach which clearly limit its applicability. 

In 0 the corrections to the BFKL kernel have been summed to obtain the corrected 
high-energy behavior of the corresponding amplitudes. This summation resulted in a 
multiple integral over intermediate rapidities, which was performed in the saddle point 
approximation and using the asymptotic form of the Is order (leading-log) BFKL Green 
function. Even after these simplihcations this leads to a highly nontrivial multiple 
integral, the calculation of which required much skill and ingenuity. A similar result 
for the high behavior has been found also in [^. 

In this note we draw attention to the fact that there exists a simple way of exactly 
solving the BFKL eigenvalue equation with the kernel being known in the two hrst 
orders. The solution turns out to exist for any real value of the complex angular 
momentum j, from —oo to -|-cxo. It means that the singularities of the BFKL amplitude 
occupy the whole real axis in the j plane, so that there is no rightmost singularity and no 
intercept (the hrst hint of a “non-Regge” behavior was noted in [Q upon the observation 
of a non-usual form of the s-dependence). All solutions to the eigenvalue equation turn 
out to have non-perturbative features. They separate into “normal” pieces, which 
can be related to the lowest order solutions and behave similarly both in q and the 
energy, and “abnormal” contributions, which are non-oscillatory and proportional to 
exp(—l/a^fe). If we neglect the abnormal pieces, we obtain an asymptotic behavior at 
high s, which coincides with [^. 

Our method easily generalizes to the case when the coupling constant is summed to 
all orders to run with the momentum scale. In the lowest order in this running coupling 
the equation still can be written, in spite of the evident difficulties at small momenta. 
Solutions of this equation also exist for any real j. This indicates that the physical 
equation should have a kernel which is different from the lowest order one, not only in 
the running of the coupling but also in its functional dependence, to give a meaningful 
result in terms of a Mellin transform. When this paper was being completed, a preprint 
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by E. M. Levin appeared , in which a very similar method has been applied in order 
to solve the 2nd order BFKL equation, with the running coupling constant summed 
to all orders and being modihed in a specihc manner in the conhnement region. The 
contents of this preprint partly overlap with our Section 4. However the central point 
of our paper - the existence of a cut in the j plane along the whole real axis - is not 
discussed in [Q. 

2 Basic equations 

We start by presenting basic formulas necessary for our derivation. They all are either 
standard or can be taken from Refs. [|l], Q. We restrict ourselves to the forward case and 
to the azimuthally symmetric wave function 'ip{q), which is known to dominate in the 
high-energy limit. The pomeron equation will be written in the form of a Schrodinger 
equation 

H,p{q)=E^P{q), (1) 

where the “energy” is related to the complex angular momentum hy E = 1— j (= —uj 
in the usual notation [Q]). The “Hamiltonian” H has been calculated up to terms of 
the second order in the coupling constant = otg'- 

H = + ( 2 ) 

Proper functions of the lowest order (BFKL) Hamiltonian Lfd) well-known: 

H^u{q) = e^yHu{q), (3) 

where 

Vv(9) = (4) 

the unperturbed energies are (in units of agNc/Tr) 

4^^ + TV'Q-w) -2i/’(l) (5) 

(e[,^) = —x(7), 7 = 1/2 -|- iz/, in the usual notation |]l|). In this equation = 

d\n.T[z)/dz, and v runs from —cxo to +cxo. The wave functions (^) are correctly nor¬ 
malized: 

= J dAq^l)l{q)'tpu'{(l) = ( 6 ) 
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The second order part of the Hamiltonian can be conveniently written in terms of 
its action on the proper functions of the hrst order Hamiltonian |l|. Namely 


H'lpuiq) = [ksc{j^) + Kr{u,q)]^p^{q). 


The hrst term in the bracket is scale invariant: 


(7) 


( 8 ) 


where the function c(z/) can be found in [^. The second term provides for the running 
of the coupling: 


«r(z^,g) = [-2ash ln(g/p)] 

where b = (llA^c — 2A^/)/(127r). We put /i = 1 for simplicity. 
Use is independent of q. So it simply shifts the energy 

a.Nr, 


(9) 


hi) 


i 1 


dvr 


yi')) ■ 


( 10 ) 


As a result the lowest energy level at = 0 acquires a factor 


ctsNc , s a.Nr 

1 - 0 ~ 1 - 6.562 

dvr 71 

for Nc = 3 and Nf = 2, which is so disturbing from the point of practical application 
of this formalism to present momentum scales. 

The second term Kr depends on q and evidently changes the equation itself. In ^ 
the inhuence of this factor was investigated by summing all orders in its action. We 
are going to propose a different approach. 


3 The z/-representation 

Our idea is quite trivial (and not new, see 0). Since the action of the second order 
Hamiltonian on the proper functions Tpi, is known, we are going to pass to the repre¬ 
sentation which these functions provide. To this aim we present any solution \l/(g) of 
the total Hamiltonian H as & superposition of 

^{q) = J duf{u)7p^{q). ( 11 ) 

Mathematically it is nothing but a Fourier transformation of q'^iq) with respect to 
\nq. Using (|^) we have 

f{iy) = Jd‘^qijl{q)'^{q). ( 12 ) 
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In the V representation the Schrodinger equation (|I]) reads 

jdu'H{u,u')f{u') = Ef{u), (13) 

where, of course, 

H{u, u') = = J (14) 

Using (^ we obtain 

=e^^^6{u-v'). (15) 

The scale invariant part in gives similarly 

T^') = -^^c(z/)e^^)5(z/ - u'). (16) 

The running coupling part involves Ing, which leads to a derivative of the ^-function, 

u') = -iasbe^J^}6'{u - u'). (17) 


We note that in all second order terms we can substitute by dehned by ([l0|) in 
the adopted approximation. 

Thus the Schrbdinger equation in the v representation turns out to be a differential 
equation of the hrst order in v. 


[ei./(^^)] + = Ef{v). 

This equation is solved trivially. The solution can be taken as 


C 


iE 


dv' 


IV 


= — exp — /-- , 


a.b Jo e^i 


oisb. 


(18) 


(19) 


where C is a constant. Note that the 1st order has zeros at z/ = i/q — ±0.6375. So, 
in this approximation, one has to specify the rule to circumvent the pole singularity in 
the integral, say, change ^ — iQ. At the 2nd order does not vanish for real v, 

so that this change is not necessary. However in both cases contains an imaginary 
part. Correspondingly, for the conjugate equation one has to take = e*. Then it 
will have similar solutions f ^ as (^) without the prefactor l/e^^ and with an opposite 
sign of Imci/. Imaginary parts of 6,^ and ly will produce certain real factors in Je and 
Je, of a non-perturbative nature (with l/ag in the exponent), the product of which 
gives unity. Because of that in a product /e/e; which is of practical importance, these 
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products always cancel. Due to this circumstance, in future we shall neglect the real 
factors generated, for real u, by the imaginary part of 

One can see that Je and form an orthonormal system: Indeed 


(/eI/e') 




du 

— exp 
eu 


n{E'-E)z{u) \ 

V Oisb / ’ 


where we have dehned 



and used the fact that in f ^ appears z{i') = z*^^). Evidently 


( 20 ) 


( 21 ) 


(/eI/e') 




dzexp 


/ i{E' -E)z \ 
V asb j 


2TTasb\C\^5{E - E') , 


( 22 ) 


and choosing 101 = {2T:asb)~^/‘^ we obtain a correct normalization. 

The most remarkable property of the solution is that it exists for any real value 
of the energy E. In the next Section we will show that there exists a normalizable 
eigenfunction for any —cx) < E < +cx). Some doubts may arise in view of the 

pole singularity at z/ = z/q generated by the denominator However, inspection 

shows that at this point the exponential provides a rapidly oscillating factor which en¬ 
sures the convergence of integrals involving /(z^). So the spectrum of the Hamiltonian, 
with the second order correction included, extends from — cx) to -|-cxo. Such a dra¬ 
matic change in the spectrum is due to a highly singular character of the “interaction 
term” proportional to a Inq. From the mathematical point of view it is an operator 
unbounded from below. With the rest of the Hamiltonian bounded from below, this 
shifts the spectrum to arbitrary large negative values of E. In the Ing space it is a 
linear potential which cannot be considered as small perturbation irrespective of the 


magnitude of the coupling constant which it accompanies. The solution (|T^ is ac¬ 
cordingly non-perturbative, the coupling constant appearing in the denominator of the 
exponent. There seems to be no simple one-to-one correspondence between the pertur¬ 
bative (i.e. 1st order) and non-perturbative (i.e. 2nd order) wave functions. A detailed 
discussion of the connection between perturbative and non-perturbative solutions will 
be given in the next Section. 
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4 Transition to the q space. The high-energy limit 


In the q space the found eigenfunctions are given by Eq. (0). We obtain 

C 


'^E{q) = 


71 . 


q\/2 ■ 


du , 2 -1 iEz(i/) 

— expUz/lng — ii /—- + 


a, 6 




(23) 


with z{u) dehned by (pT]). Since is supposed to be small we can study the integral 
in z/ by the saddle point method. The same method will also provide us with the 
asymptotics of at very small and very large q. The saddle point is determined by 
an equation 

Ez' — 1 + ash In (f = 0, 


which determines v for a given eigenvalue E. It follows 


E = ej^(l — CKsfelng^). 


Putting this into (p3[) one gets a crude asymptotic estimate 


(24) 




1 


■ exp 


a.sb 


higb |i^-e„z(0]|. 


(25) 


where u should be determined from ([2^). Eq. is nothing but the usual expression 
for the pomeron energy E as a function of u with a running coupling constant to the 
2nd order asiq"^) = q;s( 1 —a^folng^). However now we have to consider it as an equation 
for z/ at a given E. 

Let us discuss the behavior in q of the solutions, for a fixed value of E (see Fig. 
1). In all the arguments presented below, ck^ is small enough (smaller than ~ 0.05, in 
order to avoid the complications found in P]) for the minimal value of e^, emin, to be[| 
at iz = 0. 

First we take E > 0. For Ing^ < l/{asb) (or as{q^) > 0), (pi]) has a pair of 
solutions for real v (Fig. 2), starting from Tz/q at In q^ = —oo and going to u = ±cx) 
when Ing^ approaches l/{ash). According to (^Sf) , for such q the eigenfunctions "^siq) 
will essentially have a similar behavior in q as the unperturbed eigenfunctions fiviq), 


with the correspondence between E and u established by means of Eq. (p^l): they 
are plane waves in the variable Ing^, apart from the common factor 1/q and a certain 
distortion due to the third term in the exponent of (^^ (clearly visible in the form 
(0)). Thus we find an one-to-one correspondence between the 1st order perturbative 
^min ^ 0, if w6 ncglcct tli6 SGCond order correction in Cmin — {asNc/ir) 4 In 2. 
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and the 2nd order non-perturbative eigenfunctions for this part of the Ing^ space; we 
call these oscillatory pieces of the solutions “normal”. With agb small this behavior 
remains valid up to very large values of q limited by the restriction Ing^ < l/{asb). 

What happens if q becomes larger, so as l/{asb) < Ing^ < (1 —One 
can see that the picture changes and Eq. (^) will now give imaginary saddle points 
(Fig. 2), starting from z/ = ±i/2 for Ing^ = l/{asb) and approaching z/ = 0 as Ing^ 
comes near (1 — E/emin)/(,C(sb), so that the wave function becomes damped in q by 
some power factors (as for the case E < 0, Eq. (|?^) but with an opposite sign in the 
exponent). This piece of the eigenfunction which is governed by imaginary u values 
will be called “abnormal”. 

Now we take Ing^ > (1 — E/emin)/{cusb)- Eq. (|^) has now solutions for real u, 
starting from z/ = 0 and going to iz/g for Ing^ —>■ cx). It is easy to see the limiting 
asymptotics when Ing^ —>■ +cx). Taking into account the denominator in (pdf ) one 
then Ends an oscillating behavior^ 

1 / E \ 

^ bIq) ~ - exp (iuo In q^ — i -- In In ) , (26) 

q \ cagb J 

where c = This behavior is valid both for very small and very large values of q 
(and also for any sign of E, see below). 

We consider now E < 0 hxed in (P^). Again we find three pieces in the solutions, 
corresponding to different values of Ing^. For In < (1 - E/emin)/{asb), (H) will 
have solutions for real u, starting from iz/g at Ing^ = —oo and going to ±cx) for 
Ing^ —(1 — E/emin)/{c(sb) ■ This piece is equivalent to the last piece discussed in the 
case E > 0. 

For (1 - E/emin)/(,asb) < Ing^ < l/{asb), the corresponding values of the saddle 
point u will have a non-zero imaginary part. If we neglect the second order correction 
to Cj,, Eq. (p4|) will give a pair of pure imaginary points ±i|z/|, with |z/| going from 0 
at Ing^ = (1 — E/emin)/{<^sb) to 1/2 for Ing^ = l/{asb). To get rid of the singularities 
along the real axis it is convenient to pass to the integration variable z in (pS]). Then, 
for Ing^ close to l/{asb), in the complex z plane the integrand will have singularities 
at points where l/e^ vanishes, that is, at poles of e^. They occur at pure imaginary 
z/ = Ei/2, ±3z/2,.... Function ^ at these points will also take pure imaginary values 
fc = 1, 2,..., with < .... In the strip |Im z\ < the integrand 

^Actually a sum of contributions from the two points vq , which differ only in sign if we neglect the 
2nd order correction in should be taken. 




will be analytic, so that the integration contonr in z can be freely shifted np and down 
parallel to the real axis. The solntions of Eq. for E < 0 stay inside this strip, 
tending to its bonndaries a.s E ^ —oo and correspondingly z/ tends to ±i/2. So one can 
always shift the integration contonr to pass throngh the saddle point. Note that if one 
takes into acconnt the 2nd order correction to ej, then the hrst order poles at n = ±i/2 
are changed to third order poles at the same point. As a conseqnence, the saddle points 
will acqnire a real part and tend to ±z/2 at certain angles when E —>■ —oo, which will 
however not inflnence the hnal resnlt. If we take as very small, with q in this range, 
the two saddle points approach their limiting values ±i/2. The prodnct euz{v) will 
be also pure imaginary with the same sign and its modulus greater than |z/|, since \ey\ 


grows towards v = Then from (|^) it follows that we should shift the integration 

contour down to pass through the saddle point with a negative imaginary part. As a 
result we shall obtain 


~ — exp <! - ( ^ -\nq^] \u - euz{i' 




ash 


(27) 


which shows that at hnite q the abnormal pieces are damped by the non-perturbative 
damping factor exp(—const/cts)- This piece will correspond to the second piece dis¬ 
cussed for E > 0 (but, as mentioned above, with the opposite sign in the exponent: 
for E < 0 Ing^ approaches l/{ash) from below and for E > 0 from above, see Fig. 1). 

If we are now interested in the region Ing^ > 1 /(q;s 6), then the saddle point starts 
from ±cx) at Ing^ = l/{ash) and tends, for Ing^ —> -|-oo, to the point Euq at which 
e„ = 0. As a result one obtains the same oscillating asymptotic behavior as for the 
hrst and third pieces in E > 0. 

Excluding the region Ing^ > l/{asb), which for small values of as covers only 
exceptionally large values of q, we can summarize the situation saying that normal 
solutions (pieces) have an oscillating behavior and abnormal ones are power damped 
both for small and large values of q and also contain a non-perturbative damping factor 
exp(—const/tts). Abnormal solutions have no correspondence with the perturbative 
ones. Since there exist normalizable solutions for arbitrary large negative values of 
energy, that is, for arbitrary large positive values of the angular momentum j, they 
will lead to contributions to the amplitude which grow inhnitely fast at high energies. 

In fact even the passage to the energy representation results impossible due to their 
existence. The Green function of the total Hamiltonian as a function of the angular 
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momentum j is given by the spectral representation 

G{3,quq) = / dE— - , 28 

J 3-1+E 

where the integration runs over all the spectrum. The Green function as a function of 
rapidity Y = Ins is obtained by integrating (^) with a weight exp [Y (j — 1)], 




(29) 


Evidently this integral is ill-dehned for E > 0 if the integration goes from arbitrary 
large negative values of E, as in our case. 

However, throwing out non-perturbative contributions which behave cx exp(—l/a^), 
one obtains an apparently reasonable asymptotics, which coincides with the one found 
in 0]. Let us cutoff the integral over E in (|2^) by some negative lower limit Eq < emin, 


where e^m is the discussed minimal value of energy E for which Eq. gives real 
solutions for u: 

poo _ 

G{Y,q,,q)= dEexp{-YE)'^E{qi)^*Eiq)- (30) 

J Eq 

Putting our solutions we get then 


G{Y,qi,q) = 




2Ti‘^qqi Jeq 


dE exp{-Y E) 


dudui 


exp[zz/i Ing^ — iulnq^ — ifd{vi — u) + ijdE^Zi — z)], 


(31) 


where we denoted f3 = l/{asb) for brevity. We integrate over E to obtain 

<^(^,^ 1 ,^) = exp(-yEo) f [Y - if3{zi - z)]~^ 

27i^qqi J 

exp[ii/i Ing^ — iizlnq^ — ij3{vi — v) + if3Eo{zi — z)]. (32) 

Now we pass to the variable zi. The integrand in (p2D has an explicit pole in 
it, and also singularities along the imaginary axis due to the singularities of ui(zi). 
The latter are situated at hnite distance from the real axis on both its sides. The 
explicit pole, on the contrary, is quite close to the real axis, due to smallness of ag- It 
lies slightly below the real axis in the zi plane. With Eq < 0 we can close the contour 
around the singularities in the lower zi semiplane. The contribution from the cut along 
the negative semiaxis will contain a damping factor exp(—l/a^). We shall neglect it, 
which is certainly true in the limit > 0 and amounts to throwing out all anomalous 
solutions. Then we are left with only the residue at 


Y — iP{zi — z) = 0. 


(33) 
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In this approximation the dependence on Eq disappears and we get an integral over v. 

ICP r 

G{Y,qi,q) = - - / duexp[iuilnql — iulnq'^ — if3{ui — u)], (34) 

in which Ui should be determined as a function of u from Eq. (|3^) . Solution of this 
equation can be accomplished by perturbation theory, recalling that /3 = l/{asb) and 
is large. Then we obtain in the hrst three orders 

• 3 

where we denoted y = agbY. To further simplify we take into account that the sad¬ 
dle point u in the integration in (|3^) is small, of order 1/E. This allows to use an 
approximation 

Sy = eo + i8u + ai/'^, (36) 


Az/ = ui — u = 


• 2 // 

ly y z 


(35) 


with 


eo = 


Nra.^ 


TT 


41n2 


5 = 


1 - 


Nra^ 


25.8387 + 0.1869 ^ + 3.8442 


Nf 


47r 


f NcU, 
V TT 


Nr. 


Nf 

15.4262 - 2.8048 

Nc 
N 


a = a 


/Nras\^ f Nf 

(322.188- 3.10189 ^ + 21.6732 


/ 
m 


a = 


iV3 

c , 


Ncas 

71 


14C(3) 


((/(z) being Riemann’s zeta function), and to restrict ourselves to terms up to the 
second order in u in the exponent in (|34D P|. The exponent then becomes a polynomial 

(37) 


P(i') ^po-ipii' - P2i'^, 


where from (|55D and (|36|) we hnd (up to terms cx ag~^‘^Y^) 

Po = -Yeo{l - asblnql)+ ^{asbeofaY^-^dasbeoY'^, 
q^ 

Pi = In ^ — agbaeoY'^ + 6Y, 

Ti 

P2 = oY. (38) 

Integration over u gives the desired asymptotics at large Y: 

_“’‘p (p" ■ is) ■ 

^The values of as considered here are those small enough to keep a > 0, i.e. as smaller than 
^ 0.05. The case a < 0 has been examined in [^, ^ and leads to an oscillatory behavior of the Green 
function. 
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Putting expressions (|38|) into (^9]), retaining only terms up to in the exponent and 
neglecting terms of order (5/a) we obtain 


GiY,qi,q) 


TiqqiP V aY 

1 - ^(lng2 + lng2) 



(40) 


ln^(9VgD 

4aY 




This expression coincides with the asymptotics found in [^], if we take into account 
that our a = AD in [Q. 


5 Running coupling in all orders 


An amusing application of the described formalism is a possibility to write an infrared 
stable equation for the case when the interaction is taken only to the hrst order in the 
coupling, but the coupling is taken running in all orders. In other words 


where 


H'tpv = 

(41) 

eWfol = 

TTOln ^ 

(42) 


and is dehned by (|). In the following we put A = 1. 

Of course, we understand that (^2]) is quite unsatisfactory from the physical point 
of view: it implies that the coupling constant is prolonged as a pure imaginary quantity 
into the conhnement region g < 1. So the following serves only as an illustration of 
the power of the employed technique and also as a warning against a simple-minded 
use of the saddle point approximations. Note that in [^, where, as mentioned in the 
Introduction, a similar problem is solved by a similar technique, a different choice was 
made: 


Ing^ 


r(g). 


where r{q) was to be determined from a transcendental equation 

r = Ing^ - ln(6r). 

However for < \j2e/h this equation leads to complex values of r and consequently 
for as{q^) = l/(6r), which seems to be even worse than our choice. 
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In the u representation we obtain the Hamiltonian as 


Nr. 


v') = ^ie[,!^sign(i/' - u). 


The eigenvalue equation becomes 

N, 


2'iib 


Differentiating with respect to u we obtain 


Nr 


i j - r)f(v') = 




with a solution 


ei ^ C' 
.fit') = gCXp 


N,. 


dlde-l'' 


( 43 ) 


(44) 


(45) 


(46) 


TibE Jo 

This seems to be a valid solution for any real positive or negative. The conjugate 
equation is 




Tib 


(47) 


with solutions similar to (^) but with an extra factor They also exists for any E. 
Requiring (/£;|/e') = S{E — E'), we get \C\ = {2tt%/Nc)~^^‘^. Solutions (H) and its 
conjugate can be readily expressed in terms of the T-function. Indeed 

r(l/2 + iv) 


^(^) = / dv'=—i\n ■ 
Jo 


— 2iy 


r(l /2 — iu) 

In the momentum representation the eigenfunctions have the form 

C f ( N \ 

^ / dz^exp ( iulnq^ — i —) . 
nEax/2J ^ 2nbE) 


(48) 


(49) 


7rEqy/2 ■ 

To find the asymptotics at large |lng^| we employ the saddle point approximation. 
The saddle point is determined by a simple equation 

Nr 


Ing" = 


-e(') 

nbE '' 


or 


E = 


Nr 


41) 


where we used the definition (^^. This equation is again nothing but the relation 
between the energy and v, expected from m- It has to be considered as an equation 
for the saddle point u for a given E. We will now analyze, as in the previous Section, 
the behavior of the solutions in momentum space, for fixed E (having in mind that the 
saddle point will only give a reliable solution for very large | lng^|). 
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Let us take E > Then, for Ing^ < {'nhE) = —4iVcln2/(7r6i?), (|50|) 

has solutions for imaginary z/, going from v = ±i/2 at Ing^ = —cx) to z/ = 0 at 
Ing^ = Nce^^^/{7rbE). As a result, for we obtain a falling asymptotics (see below 
the case < 0); this would correspond, in the language of the previous Section, to an 
abnormal piece. 

Now, if Ing^ > Nce^i^/i'^hE), (^OD has solutions for real z/, going from z/ = 0 at 
Ing^ = to z^ = ±cxd for Ing^ = +cx). This will result in an oscillating 

asymptotics for (the normal piece), apart from the dimensionful factor 1/q. 

We consider now the case E < 0. For Ing^ < {irhE)^ the saddle point 

equation (|^ will give real solutions, going from z/ = ±oo for Ing^ = —cx) to z/ = 0 
at Ing^ = {nbE). The solution (^^ will be, as for the second piece in the case 

E>0, a normal oscillating piece. 

For In > N/^Ji^bE), (13) will have a pair of conjugate purely imaginary 
solutions for v. At large Ing^ they approach points -Ei/2. In the vicinity of, say, 
u = i/2 we have 


di) 


~ —(1/2 + iv) 


-1 


so that from (^0]) we hnd 


+ ii> = 


Nr. 


(51) 


2 ' " TibElnq^ 

Function has a logarithmic singularity at this point so that it will be approximately 
given by 


z ~ i In 


If 


= z In ■ 


Nr 


(52) 


'Kb\E\ In q^ 

where we have taken into account that in the considered region < 0. It follows that 


at this saddle point the exponent in (|49|) becomes 


zz/lng 


Nr 


2 


TibE 


~ — Ing^ 


Nr 


Eb\E\ Ing^ 


1 — In 


Nr 


'Kb\E\ In q^ 


At the complex conjugate saddle point, in the vicinity of zz = — (l/2)z, the exponent will 
have an opposite sign. In the complex v plane the integrand in (|49|) has singularities 
along the imaginary axis for |Im v\ > 1/2. In the strip |Im v\ < 1/2 it is analytic, so 
that the integration contour can be freely shifted above or below the real axis in this 
interval. Depending on the sign of 


1 1^1 ] ^ 

2 'Kb\E\\nq^ y 7r6|i?| Ing^y ’ 
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one can always shift the contour to pass over one of the two conjugate saddle points so 
as to have a negative coefficient before Ing^. As a result we obtain a falling asymptotic 
(abnormal) behavior in this Ing^ range (corresponding to the first piece in the case 
E > 0): 
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(53) 


Summarizing, both for > 0 or A' < 0 we find that the solution has two pieces, 
normal and abnormal, depending on the In range we are studying. The eigenfunction 
(^) is normalizable for any real value of E. Restricting ourselves to the case In > 
(corresponding to Ing^ > 0 for large negative E), we hnd a normal 
oscillating solution for > 0 and an abnormal falling solution for < 0. 

With the spectrum of the Hamiltonian extending from —oo to oo, we again 
have grave problems in passing from the angular momentum to rapidities. From the 
start the integral (^), which determines the Green function as a function of rapidity, 
is badly divergent at large negative E. It is amusing that nevertheless one obtains a 
seemingly reasonable asymptotics at high Y and momenta if one forgets about the ini¬ 
tial divergence and makes some crude approximations. Indeed, with the eigenfunctions 
given by (|^), we have 


G(Y,qi,q) 


exp (iui Ing^ — iulnq^ — — z)) , 

\ ■noE J 


(54) 


where Zi = zivi). The integral (|5^) does not exist. However let us forget this and 
calculate its asymptotics for large Ing^ and \nq\. Then small values of v and vi and 
consequently of z and zi evidently give the dominant contribution. At small v 


z = 





so that we obtain 


G{Y,qi,q) = 


2|Gp6, 


2dl) 


^ exp(-F£:) (5 I 
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This integral exists, since the ^-functions ensure that the integrand is zero nearly for 
all E and for large negative E in particular. We obtain in this manner 


^ \G\^Tib S{q-qi) 

G{Y, gi, q) = —— -exp 

Nr. q 


-Y 


nb In 


(55) 
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The behavior in Y is just what one would expect from naive expectations of having a 
running coupling in the intercept. 

Of course, this exercise has very little meaning. It only shows that using poorly 
controlled approximations one can arrive at seemingly sensible results even though the 
exact expression has no meaning at all. In our case it shows that the gluon interaction 
should have a different functional dependence to cure its bad qualities, which arise 
when one simply changes a hxed coupling by a running one. In mathematical terms 
the interaction then results not bounded from below, which shifts the spectrum down 
to —oo. 


6 Conclusions 


Different aspects of the 2nd order corrections to the BFKL pomeron [Q, ^ have been 
studied recently [^, |^, |] using different techniques. Some of the resulting features, as 

the existence of a “non-Regge” term oc exp (Y^) |^, or of oscillatory solutions p, ^ 
for Os greater than ~ 0.05, together with the diffusion into the infrared region 110. 
impose severe limits on the possible applicability of these results to the experimental 
situation. 

In this note we hnd the solutions to the 2nd order BFKL Hamiltonian and analyze 
the resulting spectrum in the energy plane {E = 1 — j = —u). It turns out not to 
be bounded from below, and we show that for any —cx) < E < +cx) there exists a 
normalizable eigenfunction This solution exhibits non-perturbative features. 


with the coupling constant appearing in the denominator of the exponent, Eq. (pS]). 
Both for the cases E > 0 and E < 0, solutions to the eigenvalue equation have three 
different pieces in momentum space, depending on the relation of the value of Ing^ 
with as and E. The different pieces can be classified as being “normal”, with an 
oscillating behavior as in the 1st order case, and “abnormal”, which are non-oscillatory 
in q and have no equivalence to the 1st order solutions. For very small the case 
E > 0 is almost entirely normal, whereas for E < 0 the eigenfunction is almost 
everywhere abnormal. From our study we conclude that, in the presence of the 2nd 
order corrections, the unbounded spectrum of the kernel makes the use of the complex 
angular momentum formalism to study the high-energy behavior somewhat doubtful. 

A calculation of the high energy behavior of the 2nd order BFKL Green function 
has been presented in [Q. In this paper, Mellin transforms were applied only to the 
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leading-log BFKL Green function, and the problem of the unbounded spectrum of the 
second order corrections was avoided. So the iteration techniques used in [|] seem to 
offer a possibility to compute the high energy behavior of the 2nd order BFKL Green 
function bypassing the difficulties of the unbounded spectrum found in our paper. We 
have shown that with a modification of the Mellin transform the Green function can be 
defined by simply excluding the abnormal pieces of the eigenfunctions; this definition 
leads to a high energy behavior which agrees with |^. 

In also the high-energy behavior of the 2nd order BFKL Green function is 
investigated. The result is consistent with Ref. Q]. Although the methods used in ^ 
are similar to parts of our study, no particular attention has been paid to the energy 
spectrum of the 2nd order BFKL kernel which constitutes the main purpose of our 
paper. 

So far we have concentrated on the spectrum of the 2nd order BFKL kernel, as 
it was presented in |Q. In the context of the leading-order BFKL equation |^, ^ 
it has already been emphasized that the physical spectrum of the Pomeron depends 
upon the (nonperturbative) infrared behaviour, i.e. upon the (unknown) extrapolation 
of the perturbative BFKL kernel into the small-g region. We have not yet adressed 
this question, but we feel that the arguments which in p, ^ have been given for the 
leading-order BFKL equation should now be reconsidered in presence of the 2nd order 
corrections which lead to this dramatic change of the energy spectrum. 
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Figure 1: Behavior in q of the solutions (^) for hxed E > 0 and E < 0. Real 
V correspond to an oscillatory (normal) piece, while imaginary v correspond to an 
exponentially damped (abnormal) one. 
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